Upper limits to surface force disturbances on LISA proof-masses 
and the possibility of observing galactic binaries 
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We report on the measurement of parasitic surface force noise on a hollow replica of a LISA 
(Laser Interferometer Space Antenna for the observation of gravitational waves) proof-mass (PM) 
surrounded by a faithful representation of its in-flight surroundings, namely the capacitive sensor 
used to detect PM motion. Parasitic forces are detected through the corresponding torque exerted 
on the PM and measured with a torsion pendulum in the frequency range 0.1 30 mHz. The sensor 
electrodes, electrode housing and associated readout electronics have the same nominal design as 
for the flight hardware, including 4 mm gaps around the PM along the sensitive laser interferometry 
axis. We show that the measured upper limit for surface forces would allow detection of a number 
of galactic binaries signals with signal to noise ratio up to « 40 for 1 year integration. We also 
discuss how the flight test under development, LISA Pathfinder, will substantially improve this 
limit, approaching the performance required for LISA. 



Introduction. LISA is a space-borne gravitational 
wave detector under development by the American 
(NASA) and European (ESA) space agencies [1]. The 
detector measures, by means of a laser interferometer, 
the gravitational wave-induced change in the distance 
between free-falling proof-masses (PMs) separated by 
R:! 5 X lO^m. Three PM pairs define three lines that form a 
nearly equilateral triangle. The PM pairs are assumed to 
be in pure geodesic motion, free of spurious relative accel- 
eration noise along the measurement axis (referred to as 
X from now on) to within \/2x {3x 10~^^(m/s^\/Hz)(l -|- 
(//Smnz)**)^/^}, where / is the frequency. 

The PMs in LISA are cubes made of a Au-coated Au- 
Pt alloy and as such are solid conducting bodies with 
no mechanical contact to the surroundings. The relevant 
disturbing forces can be divided into forces into those 
that act on the surface, which are most significant, and 
those acting on the bulk, namely low frequency magnetic 
and gravitational effects. Both the magnetic fields gen- 
erated inside the spacecraft (SC) and the PM magnetic 
properties can be measured on ground with sufficient ac- 
curacy, and thus the resulting magnetic force noise ex- 
pected in flight can be well characterized. The interplan- 
etary magnetic field is known to be stable enough to be 
a minor entry in the LISA error budget [2]. As there 
are no continuously moving masses in the current LISA 
design, the most worrisome gravitational noise source is 
the SC thermo-mechanical distortion. Were this distor- 
tion known in detail, the gravitational noise on the PM 
would then be well modeled based on the measured mass 
properties. Modeling of thermo-mechanical distortion is 
part of the LISA Pathfinder mission [3], the LISA pre- 
cursor under development by ESA, as is an in-flight mea- 
surement. 

Surface forces include, among other disturbances, low 
frequency electric field fluctuations, back-action from the 
electrostatic position readout, thermal radiation effects, 
and residual pressure gradients. These forces are the 



most difficult to model, as they may be connected, for 
instance, to electrochemical or radiative properties of the 
PM and electrode surfaces. 

In LISA the PM is surrounded by a system of elec- 
trodes used to detect the PM motion within the SC. 6 
pairs of electrodes are used to sense the proof mass mo- 
tion in all degrees of freedom, using 100 kHz capacitance 
bridge readouts. In addition a set of biasing electrodes 
on opposing faces along one or both directions normal 
to X provides the PM ac-voltage polarization required 
by the bridge to work. The position information is fed 
back to control loops that suppress relative SC-PM mo- 
tion. We call this electrode system, together with its sup- 
port housing, the necessary readout coaxial cables, and 
other accessories, the Gravity Reference Sensor (GRS). 
The GRS creates a nearly closed cavity around the PM, 
and, with its Au-coated electrodes and housing, serves 
as an electrostatic shield. Thus the GRS shields the PM 
from external sources of surface force and is, in practice, 
the dominant source of surface force disturbances. 

We proposed [4] a torsion pendulum to investigate sur- 
face forces on LISA-like PMs. A flrst measurement cam- 
paign [5] employed a prototype sensor with a hollow, 
40 mm Au-coated Ti PM, suspended in-axis with a W 
wire, with 2 mm gaps to the surrounding sensor electrode 
surfaces. The hollow PM allows use of a thinner torsion 
fiber, and thus higher force sensitivity, than with the full 
2 kg LISA PM, without sacrificing representativity for 
surface forces [4, 5]. The electrodes were Au-coated Mo 
plates, with biasing electrodes only along one axis (the 
z-axis). A simple home-made sensor readout electronics 
was used for the measurements. 

In connection with the implementation of the LISA 
Pathfinder mission [3] , the design of the GRS has evolved 
[6]. The PM is now a 46 mm cube, separated from the 
electrodes by a gap of 4 mm along the a:-direction, and 
2.9 mm and 3.5 mm on the other axes. The electrodes 
are now Au-coated ceramic pieces mounted inside a Mo 
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FIG. 1: Left: the electrode housing of the GRS used for the re- 
ported measurements. The AU-coated ceramics electrodes for 
the x-direction motion readout are the pair of shiny vertical 
plates visible inside the sensor. Also visible, on the adjacent y 
face, is a y-sensing electrode, the upper horizontal plate, and 
one of the electrodes used to polarize the PM for capacitive 
readout (with a circular hole in the center). Right: the hollow 
PM used for the torque measurements. 



structure. In addition the biasing electrodes are present 
on all faces except the ones along x. 

In this article we describe the results of a measurement 
campaign performed with the torsion pendulum on a pro- 
totype of this new sensor design, connected to a read- 
out electronics closely following the final fiight-hardware 
design. From the results we derive an upper limit for 
surface forces on the PM for LISA. Finally, we briefly 
compare this performance with that expected from LISA 
Pathfinder for the entire force disturbance budget. 

The main features of the prototype are described 
above, with a photograph shown in Fig.l. Details of the 
readout scheme can be found both in [5, 6]. Also shown 
in Fig.l is the pendulum inertial element, featuring a 
46 mm hollow Al cube suspended by a W wire by means 
of an Al shaft. The shaft is electrically isolated from the 
PM via a quartz ring hidden underneath the top face of 
the PM. A mirror used for the independent optical read- 
out of the angular motion is attached to the shaft, along 
with a structure used to balance the entire setup. The 
total moment of inertia lo = (4.31 ± 0.01) x 10~^ kg m^. 
The PM is suspended under high vacuum (10~^ Pa). The 
resulting torsional oscillator has a period approximately 
To = 564 s, the exact value depending on the applied 
electrical field, and a quality factor Q = 2900 ± 500. 

The GRS is mounted on a micro-manipulator that al- 
lows moving and centering it around the PM, to which it 
has no mechanical contact. This manipulator, together 
with an additional one at the upper pendulum suspen- 
sion point, allows the centering in all degrees of freedom 
with fjbm precision. The pendulum motion is measured 
both by the GRS and by an independent optical read- 
out based on a commercial autocollimator. The GRS 
has a sensitivity of order 2 nm/\/Hz in displacement and 
200 nrad/-\/Hz in rotation in rotation down to frequen- 
cies of 1 mHz, with no reliable measurement below this 
frequency. It provides the PM motion readout in all 6 
degrees of freedom. The autocollimator has an intrin- 
sic resolution of 20 nrad/VHz within the same frequency 



range and provides the readout for both the torsion angle 
and pendulum tilt around one horizontal axis. The GRS 
is equipped with heaters and thermometers, while optical 
fibers can shine UV light on PM or electrode housing to 
discharge the PM. A 3-axis magnetometer monitors the 
low frequency magnetic field fiuctuations with a resolu- 
tion of 1 pT/\/Hz above 1 mHz and a 1// upturn below 
this frequency. More details on the pendulum can be 
found in [4, 5]. 

The angular motion of the pendulum (f){t) is converted 
into an instantaneous external applied torque N{t) as 



N{t)=Io{'Ht) 
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In practice, the derivatives are estimated from a sliding 
second order fit to 5 adjoining data. Eq.(l) approximates 
the damping as viscous, while in reality it is dominated 
by structural dissipation within the fiber. However, the 
difference between these two models has no detectable 
effect on data processing. 

In order to get an estimate of the expected read- 
out background noise, we formed the average, N{t) = 
WcRsit) + Nac{t)] /2, and semi-difference, 5N{t) = 
[NGRsit) - Nac{t)]/2 of the torques, NcRsit) and 
Nac{t), calculated, respectively, from the GRS and au- 
tocollimator angular readouts cpoRS (t) and (l)ac {t)- For 
ideally uncorrelated readout noise, the power spectral 
densities (PSD) of these signals should be: 
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where Sx if) is the PSD of process X at frequency /, 
Sn.4,cns if) ^-iid ^n,4>aa if) ^^'^ the PSDs of the uncorre- 
lated angular readout noise of, respectively, the GRS and 
autocollimator, and H {f) = {loi^irf [{l/Tof - f^ + 
i if /ToQ)]}^^ is the pendulum torque-to-angle transfer 
function. Thus within this model one can estimate the 
torque PSD by: 

Sn if) = % (/) - SsN if) = Re {SN^,,Nans if)} (4) 

Eq.(4) indicates that in this model the difference 
Sn if) ~ SsN if), attributed to true torque noise mea- 
sured above the readout limit, is equivalent to the cor- 
related noise between the two readout signals, and can 
be calculated as the cross-power spectral density (CSD) 
R-ejSjv^e.ATGHs (/)} between Nqrs it) and Nac it). Note 
that, in the absence of a definite excess torque, the result 
of this estimate, if calculated as a PSD difference or GSD, 
can give negative values due to statistical fluctuations. 

PSD and cross-spectral densities are estimated by stan- 
dard Welch periodogram, using a Blackman-Harris 3'''^ 
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TABLE I; Results of the estimate of excess noise coefBcients. 



f [Hz] 

FIG. 2: Total torque noise. Data points: square root of torque 
noise PSD estimated from autocoUimator and capacitive sen- 
sor cross-correlation after time domain subtraction of the ef- 
fect of PM motion along y. Values and error bars are esti- 
mated as the mean and its standard deviation of PSD data 
falling within the corresponding frequency bins. Bins have 
equal logarithmic width. Dashed lines delimit the error band 
for the estimate of thermal noise. The continuous smooth line 
is a linear least square fit to data as described in the text. 



order window and by averaging over 20 data segments, 
adjoining data segments overlapping by 50%. Data seg- 
ments are individually de-trended by subtracting the best 
fit straight line to the data. The accuracy of the method 
has been tested with a variety of simulated data. Sim- 
ulation shows that in the presence of large drift or in 
general substantial noise power below the measurement 
bandwidth, the two lowest frequency points follow the 
frequency shape of the side lobes of the spectral window. 
These points are then discarded during further data pro- 
cessing. In Fig. 2 we report an example of the outcome of 
this data reduction procedure. 

In this specific run, following an air re-entry into 
the vacuum chamber, an extra coupling to motion of 
PM in the y-direction was detected. N (t) data conse- 
quently showed a very strong correlation with this mo- 
tion as measured by the GRS. This motion is a few 
orders of magnitudes larger than any SC motion rele- 
vant for LISA, and thus a small coupling may convert 
into a large torque that would not translate into a cor- 
responding force on a LISA PM. We thus subtracted 
the 2/-signal of the GRS multiplied just by a constant 
from the N (t) data series. The constant was chosen 
by minimizing the residual noise at low frequency and 
was found to be dV/dy = (-9.8 ± 2.2) nN. This num- 
ber, confirmed in other runs, agrees with the results of 
coarse measurement obtained by moving the GRS with 
the micromanipulator that gives dV/dy ~ 7 nN. Coin- 
cident with the appearance of this extra coupling, addi- 
tional angular stiffness was detected, measured by rotat- 
ing the GRS [8]. This extra angular stiffness amounts 
to d'^T/dclP = (1.43 ± 0.02) x 10"!° Nm/rad. The data 
for dr/dy and for d^T/dcj)^ taken together, suggest a lo- 
calized charged feature on the j/-face at ~ 1.5 cm from 
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f-^coeff. [N^m^Hz] 


f*coeff. [N^m^Hz^^] 


46 mm PM 
4 mm gap 


(3.7 ±0.3) X 10-^'' 


(3.8 ±0.4) X 10"^^ 


40 mm PM 
2 mm gap 


(5.0 ±0.4) X 10"^'' 


(1.3 ±0.2) X 10"^^ 



the center, compatible with the observed slow decay of 
the effect with time. We suspect that such a charge is 
due to one of many observed defects in the electrode Au- 
coating that expose the potentially-charged dielectric un- 
derneath. The overall effect of subtraction is to decrease 
the noise by a factor 2 between 0.2 mHz and « 2 mHz, 
with negligible effect at higher frequencies. 

The figure also shows the expected contribution from 
the thermal torque noise associated with pendulum 
damping and providing the ultimate torque noise floor. 
As damping at these low frequencies is dominated by 
dissipation within the fiber, thermal noise is expected to 
have a PSD 



Sth if) = 4fcBr 



lo (2^/ro) 
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with T the temperature and ks the Boltzmann constant. 

The figure shows that even after subtraction of the ef- 
fect of motion along y, residual noise coincides with the 
expected thermal background, within the errors, only 
near 1-2 mHz. In order to quantify the extra noise 
at other frequencies we performed a fit to the excess 
with a simple polynomial law of frequency in the range 
0.2±30 mHz by a weighted least square procedure. Many 
different power laws were tested with linear combina- 
tions of both positives and negative exponents. The only 
choice giving fully significant F-test value for all terms in 
the polynomial is the superposition of a 1//^ law with 
a term proportional to This is true for all other ex- 
perimental runs we have investigated (see below). The 
results of the fitting are reported in Fig. 2. Coefficients 
and uncertainties are reported in Tab. (I). The same table 
reports the results for a run performed with the 40 mm 
PM, 2 mm gaps prototype. Both runs were performed 
with same W torsion fiber. 

Discussion. Fig. 2 shows that excess noise only be- 
comes negligible around 2-3 mHz. Below this frequency 
the 1//^ term is a significant contribution to the noise 
that brings total noise to a value which is 40% above the 
background at 1 mHz and 2 times this one at 0.2 mHz. 
We have not yet understood the origin of this excess. 
There is no significant correlation with the magnetometer 
readout. A weak, barely significant correlation between 
N (t) and the thermometer monitoring the vacuum col- 
umn enclosing the torsion fiber was found. This temper- 
ature was also correlated with the vertical motion z{t) 
of the PM that consequently showed a correlation with 
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N (t). A subtraction of this effect would lower the excess 
by only 10% . 

A few considerations seem to indicate that the excess is 
a property of the pendulum and is not connected to the 
interaction between the PM and GRS. First the noise 
level has not been significantly affected by the change 
in the design of the GRS (Tab. (I)), with the difference 
between the results obtained with the two prototypes be- 
ing well within the variability, of a factor w 2, observed 
for each prototype among different runs taken at differ- 
ent times during the many-month campaigns performed 
in each case. This is true even though changing 2 mm 
to 4 mm gaps represents a major change, significantly re- 
ducing, for instance, the coupling between charge patches 
[7] and significantly increasing the gas molecular conduc- 
tance around the PM. 

A candidate source for the measured torque excess is 
connected with the continuous unwinding of the fiber due 
to creep in the fiber material. In the runs presented, the 
fiber unwound with a rate of « 2 x 10^^° rad/s. The rate 
is reduced by temperature annealing the fiber under load. 
Unfortunately we are limited by practical reasons to an 
annealing temperature of only f« 60°C. It is interesting 
to note that a simple Poisson model where single angular 
slip events of amplitude Scj) contribute to both unwinding 
and the residual noise, yields a step size of S4> ~ 40 nrad 
and a rate of ~ 510~^ s~^. 

We searched for evidence of such glitches by running a 
Wiener-Kolmogorov filter on the data matched to step- 
shaped torque signals. The filter resolution for angular 
steps, resulted to be cr = 33 nrad, roughly the size ex- 
pected from the reasoning above. The filter did indeed 
detect several well identified large amplitude (3 — 5a) 
events, satisfying goodness-of-the-fit test criteria, and 
which accounted for 20% of the excess noise at low 
frequency. 

One can even entirely suppress the excess noise at low 
frequency by subtracting the « 200 highest events, with 
amplitude as low as ±2cr. However, while these results 
can be taken as an indication of angular glitches as the 
possible source for the detected excess noise at low fre- 
quency, they cannot be considered as an unambiguous 
proof and such events have not been subtracted from the 
data of Fig. 2. It is worth remembering that any data 
series can be approximated with arbitrary precision by a 
sequence of variable amplitude steps and thus a subtrac- 
tion procedure can possibly, and erroneously, cancel the 
noise entirely. 

The high frequency excess proportional to is 
also statistically significant. The excess is due to a 
residual correlation between the readouts at high fre- 
quency. The real part of cross-coherence r(/) = 
Re (50„.,0oHs(/)) /VS4>. if) if), is stiff 0.10 ± 0.02 
at 10 mHz, dropping to non-significant values only above 
60 mHz. 

It would be truly difficult to attribute this excess to any 



real torque on the PM. For instance, when converted to a 
voltage on one of the electrodes, this would amount to a 
PSD of 4 V/\/Hz at 300 mHz assuming an unrealistically 
large dc voltage of 2 V on the same electrode, all figures 
unlikely by orders of magnitude [8] . 

The dependence converts into a white rotational 
noise of the entire apparatus relative to the local frame 
of inertia by ^30 nrad/ v^Hz. Thus a possible explanation 
may be such a rotation due either to seismic noise or to 
local mechanical distortion. However a residual electrical 
cross-talk between GRS and autocollimator cannot be 
excluded. 

Though the detected excess is likely not due to the 
GRS itself and is not then relevant to LISA, we still in- 
clude it when the estimating the upper limit to the distur- 
bances acting on the PM. This upper limit also includes 
the fit uncertainty and by the uncertainty in the evalua- 
tion of thermal noise, conservatively estimated to be 20% 
and coming from the error in Q-factor measurement. 

Overall this upper limit for the data of Fig. 2, defining 
fo= ImHz is: 

3^,„-ig:{o(^)%0.3(ii)+1.2.I0-.(i)' 

+^10-(f)' + 2^2.,0-»(i-)'| 

These torque data can be converted into an equivalent 
acceleration noise with a few simple assumptions on the 
contributing effects [5]. For instance simulations [5] give 
that randomly distributed normal forces, like those as- 
sociated with parasitic electric fields, would exert both 
torque and force noise, with a ratio of their amplitudes 
of roughly 0.4L, with L the PM size. This holds for dis- 
turbances with spatial exponential correlation lengths up 
to L/W. The ratio would decay to i/4 for a correlation 
slightly larger than L. iy/4 is also roughly the ratio ob- 
tained for an electrical disturbance directly applied to 
the electrodes of the a;-face, as would be expected from 
electronic readout back-action. Finally forces acting tan- 
gentially, like diffuse scattering of molecules, show a ratio 
of L/2. We then convert the torque to an equivalent ac- 
celeration noise that the corresponding forces would ex- 
ert on the 1.96 kg LISA PM by using the lowest ratio 
L/i, i.e. 10.7 cm. The reader should be aware though 
that obviously our pendulum geometry would not detect 
disturbances applied uniformly and normal across the x 
face or tangentially to the y and z faces, or disturbances 
acting normally onto the center of the a;-face or in a few 
other symmetric positions. 

The results of the conversion of the data of Fig. 2 into 
an equivalent acceleration noise are reported in Fig. 3. 
Within the same figure we also report for comparison 
the acceleration level that would guarantee to LISA a 
measurement of the gravitational wave signal from a few 
binary systems in our Galaxy over one year integration 
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f [Hz] 

FIG. 3: Dashed line: upper limit to pendulum torque noise 
with GRS converted into an equivalent acceleration on one of 
LISA PMs. Dots: acceleration level required to detect with 
LISA, with signal to noise ratio of 1, the signal from the cor- 
responding galactic binary system over one year integration 
time. Dotted line: projected upper limit for non-modeled 
overall force disturbances that the LISA Pathfinder mission 
may achieve. 

time. The signal amplitude is taken from [9] and the 
LISA sensitivity is averaged for source inclination and 
polarization. 

It is important to stress again the limits of this observa- 
tion. Besides the lack of sensitivity to volume forces, all 
effects due to coupling among different degrees of free- 
dom that would take place for a fully free-falling PM 
are absent here. Additionally the environment on board 
LISA will be substantially different from the laboratory 
one, with higher PM charging rate due to cosmic rays 
but with lower level of temperature and magnetic field 
fluctuations. The effects of these disturbances, however 
are well understood, have been investigated separately by 
dedicated experiments [8], and found in good agreement 
with the expectations. 

In Fig. 3 we also report the projected result of a com- 
parable reduction process on the data from the LISA 
Pathfinder mission [3]. The mission will measure, by 
means of a laser interferometer, the relative acceleration 
of two LISA-like PMs in free-fall within the same SC 
along the line joining their centers. 

Mission requirements call for a demonstration of rela- 
tive acceleration of better than 3 x 10~^^(m/s^-\/Hz)(l + 
(//3mHz)^)^/^ for frequencies between 1 mHz and 
30 mHz. The relaxation of requirements relative to LISA 
is dictated by the need of keeping the test reasonably sim- 
ple, but does not affect the GRS design. The current esti- 
mate of the error budget [10] places the sensitivity some- 
how better than requirements at « 2 x 10"^'' m/s^\/Hz, 
with the largest contributions being the cross-talk with 
motion of the PMs along degrees of freedom other than 
the measurement ones, the magnetic field fluctuations, 
and the residual coupling to the SC to gravitational and 



electrical field gradients. This last coupling, which will 
be enhanced to allow measuring the effect, can also be 
balanced between PMs in order to suppress the infiu- 
ence of the SC coupling on their relative motion. PM 
motion along all degrees of freedom is measured with 
sufficient precision to allow for a subtraction similar to 
what we did with tilt in the case of the pendulum. LISA 
Pathfinder also carries magnetometers so that magnetic 
field effect correction should be possible. If these largest 
entries are suppressed the way we have described here, 
LISA Pathfinder may be able to put an upper limit on 
non-modeled disturbances to on a single PM accelera- 
tion of w 5 X 10""'^^ m/s^-\/Hz at 1 mHz. Above this 
frequency the laser readout noise is expected to limit the 
sensitivity, while below 1 mHz noise the degradation is 
expected to be no worse than a 1//^ dependence. This 
is the rationale behind the curve shown in Fig. 3. 

In conclusion current estimates of parasitic forces 
would indicate that, even at the present level of experi- 
mental knowledge, LISA should be able to detect bright 
gravitational wave signals. If LISA Pathfinder will meet 
its expected performance, the risk of a major unpredicted 
force disturbance will almost be entirely retired and LISA 
would have basically its entire planned science objectives 
guaranteed. 
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